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Energy cascade rates and Kolmogorov's constant for nonhelical steady magnetohydrodynamic 
turbulence have been calculated by solving the flux equations to the first order in perturbation. For 
zero cross helicity and space dimension d = 3, magnetic energy cascades from large length-scales to 
small length-scales (forward cascade). In addition, there are energy fluxes from large-scale magnetic 
. . . field to small-scale velocity field, large-scale velocity field to small-scale magnetic field, and large- 

, scale velocity field to large-scale magnetic field. Kolmogorov's constant for magnetohydrodynamics 

■ is approximately equal to that for fluid turbulence (~ 1.6) for Alfven ratio 0.5 < < oo. For higher 

' space-dimensions, the energy fluxes are qualitatively similar, and Kolmogorov's constant varies as 

(N ■ d}^^. For the normalized cross helicity (Jc — > 1, the cascade rates are proportional to (1 — (Tc)/(1+(Tc), 

. and the Kolmogorov's constants vary significantly with Oc- 

< 

in 
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I. INTRODUCTION 

\ There are various phenomenological, numerical, and theoretical result on magnetohydrodynamics (MHD) turbu- 
• lence. Kraichnan jlj and Irosnikov ^2] proposed a first phenomenology for homogeneous and isotropic MHD which 
• ' connects energy spectrum and cascade rates. In this phenomenology, the kinetic and magnetic energy spectrum 
.^^ . (i?"(fc) and E^{k) respectively) are proportional to k~^/'^, and the energy cascade rate of flux 11 = (^fc))^fc'^/i?o, 

■ where Bq is the mean magnetic field. In an alternate phenomenology, Marsch Matthaeus and Zhou 0|, and Zhou 
I I , and Matthaeus predicted that energy spectrum is Kolmogorov- like {E{k) cx fc-5/3) and H = {E{k)f''^k^''^. Verma 

et al. 1^ numerically calculated the energy cascade rates of Elsasser variable u ± b (u,b are velocity and magnetic 
, field fluctuations), and found their results to be consistent with Kolmogorov- like phenomenology, rather that that of 

■ Kraichnan and Irosnikov. Frick and Sokoloff 3 studied the spectra and cascade rates in a shell model of MHD turbu- 
\ lence; they found the spectral index to be close to 5/3 in absence of cross helicity and magnetic helicity. According to 

■ their study, the helicities suppress the cascade process. Recently, Muller and Biskamp |^ and Biskamp and Miiller 
' computed the spectral index and found it to be closer to 5/3. Kolmogorov-like spectrum is also supported by recent 

O theoretical results J^JiJ2, 13, 14, 15]. 

^""^ The energy cascade rates of MHD depend on cross helicity {He — u • b), magnetic helicity {Hm = a • b, where 

a is vector potential), and kinetic helicity {Hk = u • w, where lo is vorticity). Pouquet et al. applied EDQNM 
^ approximation to study energy fluxes. For nonhelical MHD, Pouquet et al. [13 argued that the ME cascade is forward, 
• ^ '\ i.e., from large-scale to small-scale. However, in the presence of helicity, they observed that the large-scale magnetic 
^ ■ energy brings to equipartition the small-scale kinetic and magnetic excitation by Alfven effect, and the "residual 
helicity" , Hk — Hm , induces growth of large-scale magnetic energy and helicity. Pouquet and Patterson ^3 studied 
this problem using direct numerical simulation and arrived at similar conclusions. In the present paper we will derive 
the energy fluxes of nonhelical MHD using field theoretic methods. The energy fluxes of helical MHD are discussed 



^ , in a companion paper, Verma [l8l |. referred to as paper II. 
■ - - ' Since there are two fields u and b in MHD, the energy can be transferred from u to u, u to b, and b to b. The 
resulting energy fluxes due to these transfers are illustrated in Fig. 1. These fluxes have been numerically calculated 
recently by Dar et al. and Ishizawa and Hattori ^| for two-dimensional (2D) MHD turbulence. The numerical 
values calculated by Dar et al. for tic ~ and rA (KE/ME) « 0.5 are listed in Table The prime conclusions of 
Dar et al. and Ishizawa and Hattori's 2D numerical study are: (1) the ME cascades from large-scales to 
small-scales (forward cascade); (2) there is a significant energy transfer from large length-scale velocity field to large 
length-scale magnetic field; this transfer could play an important role in ME enhancement; (3) there is an inverse 
cascade of KE. In another recent work, Cho and Vishniac |^ have derived some interesting scaling relationships 
between the energy transfer rates and verified them using numerical simulations. Since the direction of total energy 
cascade is the same in 2D and SD (three-dimensional) MHD turbulence, some of the conclusions drawn by Dar et 
al. 10 and Ishizawa and Hattori based on 2D simulations are expected to hold at least qualitatively in 3D MHD 
turbulence. Therefore, in this paper we compare the numerical results of Dar et al. [l9| and Ishizawa and Hattori |2"o| 
with our 3D analytic results. 

In the present paper we have carried out the energy cascade rate calculation for MHD turbulence for the inertial- 
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range wavenumbers using perturbative field-theoretic technique. Here we assume that the turbulence is homogeneous 
and isotropic to make the problem tractable. Even though the real-world turbulence does not satisfy these properties, 
at least at the large-scales, many conclusions drawn using this assumption give us important insights into the energy 
transfer mechanisms, as will be discussed in this paper and paper II. We assume that the mean magnetic field is absent; 
this assumption is to ensure that the turbulence is isotropic. Our procedure requires Fourier space integrations of 
functions involving products of energy spectrum and the Greens functions. Since there is a general agreement on 
Kolmogorov-like spectrum for MHD turbulence, we take E{k) cx k~^/^ for all the energy spectra for MHD. For the 
Greens function, we substitute the "renormalized" or "dressed" Greens function computed by Verma p^ . After 
this substitution, various energy fluxes and Kolmogorov's constant of MHD are computed. Using the steady-state 
condition, we also calculate the energy supply from the large-scale velocity field to the large-scale magnetic field. This 
result is quite robust, and is independent of the nature of large-scale forcing. In this paper we assume both kinetic 
and magnetic helicity to be absent. The energy fluxes for helical MHD are discussed in paper II. 

The parameter space of MHD is rather large because of various energy spectra. The two well known parameters 
are: the normalized cross helicity CTc, which is the ratio of twice cross helicity and energy; and the Alfven ratio ta, 
which is the ratio of kinetic energy and magnetic energy. Calculation of cascade rates for arbitrary CTc and rA is quite 
complex. In this paper we limit ourselves to two limiting cases: (1) CTc = and whole range of ta] (2) CTc — > 1 and 
rA = 1. Strictly speaking, the parameters used in our calculations are spectral ac{k) — 2Hc{k)/{E'^{k) + E^{k)) 
and r^(fc) = E'^ (k) / (k) . Since our calculation is confined to inertial-range wavenumbers where Kolmogorov's 5/3 
powerlaw is valid for all the energy spectra of MHD, both a{k) and rA{k) can be treated as constants. Note that 
these parameters may differ from the global Uc and va- We carry out our theoretical analysis in various dimensions. 
We will show that our theoretical results are in general agreement with the simulation results of Dar et al. and 
Ishizawa and Hattori j23| . 

The outline of this paper is as follows: In section 2 we calculate various cascade rates for CTc = case. The other 
extreme case CTc — > 1 is considered in section 3. Section 4 contains summary and conclusions. 



II. CASCADE RATES IN MHD TURBULENCE: ac = Q 



In this section we will analytically compute the energy cascade rates when CTc = 0. We take the following form of 
Kolmogorov's spectrum for kinetic energy (KE) and magnetic energy (ME) 

= i^"H2/3fc-5/3 (1) 
E\k) = E^(k)/rA (2) 

where is Kolmogorov's constant for MHD turbulence, and H is the total energy flux. Another Kolmogorov's 
constant K is defined for the total energy, 

Etotaiik) = E^{k) + E\k) = E{k) = XH2/3A;-5/3, (3) 

with 

K = K^l+r-/) (4) 

With this preliminaries we start our flux calculation. 
The incompressible MHD equations are 

^ + (u-V)u = -Vp + (b • V)b + i^V^u (5) 

ot 

dh 

— + {u-V)h = -(b- V)u + ?7V2b (6) 
ot 

V-u = (7) 
V-b = (8) 

where u and b are the velocity and magnetic fields respectively, p is the total pressure, and v and /i are the kinematic 
viscosity and magnetic diffusivity respectively. To compute various energy transfers among various Fourier modes we 
resort to the energy equations, which are [la, [la 
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+^-^(k'|p|q) + 5"^(k'|q|p)] (9) 

+5^^k'|p|q) + 5'"'(k'|q|p)] (10) 

The above integrals have constraints that k' + p + q = 0(k= — k'). The equal-time correlation functions used in the 
energy equations are defined using 

(u,(p,tK(q,t)) - Py(p)C""(p,i,<),5(p + q) (11) 
(6,(p,t)&,(q,t)) = P,,(p)C''^(p,<,i)^(p + q) (12) 
(u,(p,t)&j(q,t)> = Py(p)C"^(p,i,t)5(p + q) (13) 

and the energy transfer rates 5'(k'|p|q) are defined using 

5""(k'|p|q) = -3([k'.u(q)][u(k')-u(p)]), (14) 

^^"(k'lplq) = -3([k'.u(q)][b(k')-b(p)]), (15) 

5"''(k'|p|q) = 3([k' . b(q)][u(k') . b(p)]), (16) 

5''"(k'|p|q) = -5"''(p|k'|q) (17) 

Here 3 stands for the imaginary part of the argument. Note that C"'' — because CTc has been taken to be zero. 
The above equations are based on Dar et al.'s formalism, which is a generalization of those of Pouquet et al. 0| and 
Stanisic [23| and others. In Dar et aL's formalism, the terms S'(k|p|q) represents energy transfer from mode p (the 
second argument of S) to k (the first argument of S) with mode q (the third argument of S) acting as a mediator. 
Note that in the expression for S, the field variables with the first and second arguments are dotted together, while the 
field variables with the third argument is dotted with the wavevector k. Dar et alJs formulas has certain advantages 
over those of Pouquet et al. |0| and Stanisic Some of the quantities to be defined below were not accessible 

in the earlier formalism, but now they can be calculated using Dar et al.'s formulas. In addition, the flux formulas 
derived using the new scheme are relatively simpler. 

After some algebraic manipulation it can be shown that 

5""(k'|p|q) + 5™(k'|q|p) + 5™(p|k'|q) + 5™(p|q|k') + 5""(q|k'|p) + 5™(q|p|k') = (18) 

5'"'(k'|p|q) + 5^^(k'|q|p) + ^^^(p|k'|q) + 5'"'(p|q|k') + 5'"'(q|k'|p) + 5'"'(q|p|k') = (19) 
5"^(k'|p|q) + 5"^k'|q|p) + 5"''(p|k'|q) + 5''^(p|q|k') + ^"''(qlk'lp) + 5"''(q|p|k') 

+^''"(k'|p|q) + 5^"(k'|q|p) + 5^''(p|k'|q) + 5''"(p|q|k') + ^""(qlk'lp) + ^'^qlPlk') = (20) 

These are the statements of "detailed conservation of energy" in MHD triads (when = r] = 0) [23|. 

For energy flux study, we split the wavenumber space into two regions: k < ko (inside "fco-sphere") and k > ko 
(outside "/cQ-sphere" ) . This division is done for both velocity and magnetic fields. The energy transfer could take 
place from inside/outside u/b-sphere to inside/outside u/b-sphere. In terms of S, the energy flux from inside of the 
X-sphere to outside of the F-sphere is 

where X and Y stand for u or b. The energy fluxes from inside w/fe-sphere to outside u/b-sphere can be calculated 
by earlier formalism, as well as that of Dar et al. However, the fluxes from inside w-sphere to inside 6-sphere, and 
outside u-sphere to outside 6-sphere can be numerically calculated only by Dar et al.'s formalism. In this paper we 
will analytically calculate the above fluxes in the inertial range using the Kolmogorov-like energy spectrum. 
We assume that the kinetic energy is forced at small wavenumbers, and the turbulence is steady. Therefore, 

n,"< = Hg<+H?< (22) 
Input Kinetic Energy = H^< + H^< + n^< (23) 

We calculate the energy flux H^^ using the above steady-state property. Hence the energy feed into the large-scale 
magnetic field from the large-scale velocity field could be obtained theoretically irrespective of nature of large-scale 
forcing. 
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We will analytically calculate the above energy fluxes [Eq. HU] in the inertial range to the leading order in 
perturbation series. It is assumed that u(k) is gaussian to leading order. Consequently, the ensemble average of S^^ , 
(S*^"^), is zero to the zeroth order, but is nonzero to the first order. The first order terms for S^^{k\p\q) in terms of 
Feynman diagrams are shown below: 



(5™(fc'b|g)) 




+ 








(24) 



(5'"'(fc'b|g)) = - 






(25) 



(5^''(fc'b|g)) = 




\ / \ 





(26) 



{S'\k'\p\q)) = 






(27) 



In the above diagrams the solid, dashed, wiggly (photon), and curly (gluons) lines denote {uiUj) , {bibj) ,6*"", and 
G^^ respectively. In all the diagrams, the left vertex denotes fc^, while the filled circle and the empty circles of right 
vertex represent (— i/2)Pj^^ and —iP^^ respectively. Since G"'' , G''" , C"^ , and G^" are zero when dc = 0, we have 
not included the Feynman diagrams containing these terms. When we substitute (uiUj) , (bibj) using Eqs. H11I12|I . we 
obtain terms involving {p,t,t')C^ {q,t,t'). The resulting expressions for various (S^ ^ {k\p\q)') are 



{S--ik\p\q)) 



{S'^klplq}) 



dt' [Ti {k,p, q)G'''' {k,t~ t')C"" (p, t, t')C"" {q, t, t') 

+nik,p,q)G''^{p,t - t')C^^-{k,t,t')C^^-iq,t,t') 
+Ts (fc, p, g)G™ [q, t ~ i')G™ (fc, t, i')C™ (p, i, t')] 

dt'[T2{k,p, q)G''''{k, t - t')C''\p, t, t')C''H.q, t, t') 

+Trik, p, q)G''\p, t - i')C™ (fc, t, t')C''\q, t, t') 
+Tii{k,p, q)G^^iq, t - t')G"^{k, t, t')C'\p, t, t')] 

dt'[T3{k,p, q)G'\k, t - t')C^^{p, t, t')C''\q, t, t') 

-OO 

+n{k,p, q)G""ip, t - t')C''\k, t, t')C"'{q, t, t') 
+T,2{k,p, q)G'\q, t - t')C''\k, t, t')G™(p, t, t')] 



(28) 



(29) 



(30) 
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(S'^klplq)) = f dt'[Ti{k,p,q)G''\k,t-t')C^\p,t,t')C'^^{q,t,t') 

J —oc 

+Tsik,p, q)G'\p, t - t')C''{k, t, t')C^"{q, t, t') 

+Twik,p, q)G^^{q. t - t')C'"ik, t, t')C'\p, t')] (31) 

where Ti(k,p, q) are functions of wavevectors k,p, and q given in Appendix A. 

The Greens functions can be written in terms of "effective" or "renormahzed" viscosity i^^k) and resistivity ri(k) 
(see Verma for details) as 

G"''{k,t-t') = cxp(-i^(fc)fc2(t_t')) (32) 
G''\k,t-t') = cxp(-77(fc)fc2(t-t')) (33) 

The relaxation time for C'^'^{k, t, t') is assumed to be the same as that of G"", and that of C^^{k, t, t') is assumed to 
be the same as that of G'^^ . Therefore the time dependence of the unequal-time correlation functions will be 

C'^'ik^t^t') = exp(-z/(/c)fc2(<-f')) C™(fc,i,i) (34) 
C^\k,t,t') = e^]y{~i]{k)k^{t~t'))C^\k,t,t) (35) 

The above forms of Green's and correlation functions are substituted in the expression of (5-^^), and the t' integral 
is performed. Now Eqs. (|21l I28|l yield the following flux formula for n^>(fco): 



n:<(fco) 



dk /■ dp 1 



Jk>ko {2nr 7,<,„ {27rr u{k)k^ + v[p)p^ + u{q)q^ 
[T,{k,p,q)C-^p)G^^q)+T,{k,p,q)G^^k)G"-iq)+T,{k,p,q)G"^k)G"^p)] (36) 

The expressions for the other fluxes can be obtained similarly. 

The equal-time correlation functions C^^{k, t, t) and C''^{k, t, t) at the steady-state can be written in terms of one 
dimensional energy spectrum as 

C-^k, t, t) = ^-A_^fc-('^-i)i?"(fc) (37) 

C'\kXt) ^ ^^^^k-^'-'^E\k) (38) 

where Sd is the surface area of d-dimensional unit spheres. We are interested in the fluxes in the inertial range. 
Therefore, we substitute Kolmogorov's spectrum [Eqs. (|ll2fl ] for the energy spectrum. The effective viscosity and 
resistivity are proportional to /c"^/"^, i.e., 

v{k) = (if")l/2nl/3fc-4/3j,* (39) 

r^{k) = (if")i/2ni/3fc-4/3^*, (40) 

and the parameters v* and 77* were calculated in Verma 

The d-dimensional volume integral under the constraint k' + p -f- q = is given by |2^ 



dq = Sd-i j dpdq 



^(sina)''-3 (41) 



/p+q=k 

where a is angle between vectors p and q. We also nondimensionalize Eq. 13()l) by substituting [22 



ko ko 

— ; p=—y- q^—w (42) 

u u u 



which yields 



{d-irSdJo 



duln(l/t;) / dw{vwf-\siiiaf-'^F^<{v,w) 

l-v 



(43) 
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where the integrals Fy^iv, w) are 
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- + ^.(^"^/\ ^ ^,2/3) [^2(^, w){vw)-''-ir2' + triv, w)w-^-h^' + h,{v, w)v-^-ir-/] (45) 

- ^.^2/3 + + ^2/3) H(^H-'^~^^A^ + ^6(t^, H^-^'^^a' + tl2{v, w)v-''-ir-/] (46) 

= ^.^2/3 + ^1(1 + ^2/3) w){vw)-''-ir-/ + t8(z;, w)w-''-h-^^ + tio(i^, i.)^;-'-ir:,^]. (47) 

Here ti{v,w) = Ti{k,kv,kw)/k^. Note that the energy fluxes are constant consistent with the Kolmogorov's picture. 
We compute the bracketed terms (denoted by /y> ) numerically and find that all of them converge. Let us denote 
I ^Iu>+ lit + 4"> + Ibt- Using the fact that the total flux H is 

n = n^<+n^„< +n^<+n^<, (48) 

we can calculate the value of constant , which is 

j^u ^ (J) -2/3 (49) 

In addition, the energy flux ratios can be computed using Ily ^ /II = ly^ / 1. The flux ratio IIJ^^ /II is obtained using 
steady-state condition [Eq. (|22|l ]. The values of constant K can be computed using Eq. 10}. The flux ratios and 
Kolmogorov's constants for d — i and various ta are listed in Table ITll The same quantities for ta — I and various 
space dimensions are listed in Table Hill 

The following trends can be inferred by studying Table UTI We flnd that for d = 3, H^^/II starts from 1 for large ta 
and decreases nearly to zero near ta = 0.3. All other fluxes start from zero and increase up to some saturated values; 
this implies that near ta ~ 1, all the energy fluxes become significant. Clearly, the sign of H^^is positive, indicating 
that ME cascades from large length-scale to small length-scale. Under steady-state the large-scale ME is maintained 
by n^^, which is one of the most dominant transfers near ta = 0.5. The energy flux IIJ^^ entering the large scales 
magnetic energy could play an important role in the ampliflcation of magnetic energy. In paper II we will construct 
a dynamo model based on the energy fluxes. 

The Kolmogorov constant K for d — 3 is listed in Table ITll For all ta greater than 0.5, K is approximately constant 
and is close to 1.6, same as that for fluid turbulence {ta — oo). Since II cx K~^/'^^ we can conclude that the variation 
of TA (redistribution of fluid and magnetic energy, keeping the total energy fixed) does not change the total cascade 
rate. Near ta = 0.3, the constant K appears to increase, indicating a sudden drop in the cascade rate. When ta is 
decreased further, near ta = 0.25 both zy* « and r?* « [13, or X oo. This signals an absence of turbulence for 
rA near 0.25. This is consistent with the fact that MHD equations are linear in the — * (fully magnetic) limit, 
hence do not exhibit turbulence. However, it is still surprising that turbulence disappears near ta = 0.25 itself. 

The Kolmogorov's constant K com put ed above can be used to estimate the amount of turbulent heating in the 
solar wind. Verma et al. psj and Tu '26] have put constraints on the turbulent heating in the solar wind from the 
radial variation of temperature in the solar wind. Verma et al. observed that when K k, 1, all the heating in 
the solar wind for streams with CTc ~ can be accounted for by the turbulent heating. Our theoretical value for this 
constant in the absence of mean magnetic field is approximately 1.5, larger than 1. If we take K « 1.5 for solar 
wind streams with CTc ~ 0, only a fraction, possibly around half ((l/1.5)'^/2)^ of the heating will be due to turbulence. 
However, neglect of mean magnetic field, anisotropy, helicities etc. are gross assumptions, and we can only claim 
general consistency of the theoretical estimates with the observational results of the solar wind. 

We have calculated the flux ratios and the constant K for various space dimensions d > 2.2. In Verma 14] it has 
been shown that for d < 2.2, the RG fixed point is unstable, and the renormalized parameters could not be determined. 
Due to that reason we have calculated fluxes and Kolmogorov's constant for d > 2 only. For these calculations we 
take rA = 1, which is a generic case. The calculated values are shown in Table HTll It is striking that all the fluxes 
are approximately the same for large d. In addition, /II is approximately 0.5 for all dimensions greater than 4. 

We verify that ly^ for constant zy* and r]* are proportional to d~^. In Verma |Q| we find that v*,!]* oc d~^/'^. 
Therefore, K cx d~^/^ . This result is a generalization of theoretical analysis of Fournier et al. 0| for fluid turbulence. 

In this section we calculated the cascade rates for Uc = 0. In the next section we take the other limit ac I. 
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III. CASCADE RATES IN MHD TURBULENCE: 1 



In this section we will describe the calculation of the energy cascade rates for the large normalized cross-helicity 
(ctc — > 1), and show that the cascade rates crucially depend on cross helicity. For cases with ^ 1, it is best to 
work with Elsasser variables = u± b. For the following discussion we will denote the ratio / ^jz+p) by r. 

Clearly r ^ 1. Here we limit ourselves to r^i = 1. 

The incompressible MHD equations in terms of are 

— + (z=F.V)z± -Vp + i^+V^z^ +:/_V2zT (50) 
V-z=^ = (51) 

where p is the total pressure, and v± = {v± rf) /2. Numerical simulations of Verma et al. and Dar , solar wind 
observations of Matthaeus and Goldstein , and Marsch and Tu 30] , and theoretical calculations of Verma 0, ll^ 
show that Kolmogorov-like energy spectrum is valid even for nonzero cross helicity, i.e., 

EHk) = K^^^^^.k-^'' (52) 

where are Kolmogorov's constants for MHD. The above equation was first derived by Marsch 0|. 
The corresponding equations for the energy evolution are 

^ ■ 2i.+fc2) C±±(k,i,t) + 2i._fc2C±T(k,t,i) ^ ^ 



df ^ J - (d-l)(27r)'iA'+p+q=o(27r)'^ 

[5±±(k'|p|q) + 5±±(k'|q|p)] (53) 

where 

5±±(k'|p|q) = -3([k' . zT(q)][z±(k') . z±(p)]) (54) 
and the equal-time correlations functions C^* and are defined using 

(zf(p,i)2f(q,i)) - P.,(p)C±±(p,t,i)^(p + q) (55) 



3 

{zt{p,t)zj{<i,t)) = E,,{p)C^^{p,t,t)S{p + <i) (56) 



From Eq. (|54() it is evident that in the nonlinear transfers, the modes z+ transfer energy only to z+ while z acts as 
a mediator. Similarly transfers energy only to z~ with z'^ acting as a mediator. It can be easily shown that 

5±(k'|p|q) + 5±(k'|q|p) + 5±(p|k'|q) + 5±(p|q|k') + 5±(q|k'|p) + 5±(q|p|k') = (57) 

These equations correspond to the "detailed conservation of energy" in MHD triads. 

In terms of z^ variables, there are only two types of fluxes H^, one for the z'^ cascade and the other for z^ cascade. 
In terms of S, these energy fluxes H^ are 

As described in the last section, the above fluxes are calculated to the leading order in perturbation series. To the 
first order, (S'='=='=(k'|p|q)) are 

<J —oo 

+Tu{k,p,q)G^^{k,t~t')C^^{p,t,t')C^^{q,t,t') 
+T,,ik,p, g)G±T(fc, t - t')C^^{p, t, t')C^^{q, t, t') 
+T,e{k,p,q)G^^{k,t-t')C^^ip,t,t')C^^{q,t,t') 
+Ti7{k,p,q)G^^{p,t-t')C^^{k,t,t')C^^{q,t,t') 
+Tis{k,P, q)G^^{p, t ~ t')C^^{k, t, t')C^^{q, t, t') 
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+ri9(fc,p, g)G±^(p, t - t')C^^{k, t, t')C^^{q, t, t') 
+T2o(fc, p, q)G^^(p, t - t')C^^ {k, t, t')C^^{q, t, t') 
+T2i(fc, p, q)G^^ {q, t - t')C^^ {k, t, t')C^^ (p, t, t') 
+T22(fc, p, q)G^^ {q, t - t')C^^ {k, t, t')C^^ [p, t, t') 
+T23ik,p, q)G^^{q, t - t')C^^{k, t, t')C±T(p, t, t') 

+T24k,p, q)G^^iq, t - t')C^^ (k, t, t')C^^ {p, t, t')] (59) 



where Ti{k,p, q) are given in Appendix A. 

Now we use the approximation that r is small. In terms of renormalized v matrix 



m = ( :5, ^ h (60) 



the Green's function G{k, t — t') = exp — [i>fc^(t — t')] to leading order in r is 

Gik,t~t')^ , , ' \^ ,^ \^ ^)}^' ^ ,^ >\. (61) 

For derivation and further details on the renormalized v, refer to Verma The correlation matrix C{k,t — t') is 
given by 



( C++{k,t,t') C+-{k,t,t')\ _ A., , ,,^ /^C'++(fc) C+-{k) 

\C-+ik,t,t') C—{k,t,t') ) " ^^"'^ ^>\C-+{k) C—{k) 



(62) 



The quantities C^*(fc) can be written in terms of E^{k) as 



We take C^^(fc) — (or — 1) for simplifying the calculation. We take Kolmogorov's spectrum for E^{k) [see 
Eq. inU], and 

The renormalized parameters C* , a* , ■(/i* , and (i* have been calculated in Verma Finally, the matrices G(fc, t — t') 
and G{k,t,t') can be written in in terms of renormalized parameters and Kolmogorov's spectrum. 

Now we substitute G{k,t — t') and C{k,t,t') in Eq. 1)59(1 . and keep terms only to the leading order in r. We find 
that the terms (1,4,5,8,9,10) of Eq. (|59|l vanish. The final equation for the fluxes 11* to the leading order in r are 



n 

where the integrand are 



n± = ^i£iX(i^+)3/2 



"^^"^'^ I di;ln(l/t;) / dw{vwf (sin af-^F^{v,w) 



(65) 



1 . , ^_^_9/■^a* (1 1 



+t,,{v,w)w-'-'/'-^+t,,iv,w)w-''-'/'- ^ ^ ^ 



l3*Uj2/3 ■ , / |^^*(-„2/3 ^_y^2/3) f3*wV3 



(3* \^*(t;2/3 + ^2/3) /3*u;2/3 

F . ^i3(.,H(-)--^/^^,(^+^.(.,^)^--^/^^,^^ (67) 
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We denote the bracketed term of Eq. H65|l by and compute them numerically. We find that the integrals are finite 
for d = 2 and 3. Also note that are independent of r. We calculate the constant of Eq. H65|) in terms of the 
constants are listed in Table Hvl for various values of r in d = 2 and 3. The constants depend very sensitively 
on r. Also, there is a change of behaviour near r = (7^//+)^ = Vc,; < for r < Tc, whereas inequality reverses 
for r beyond Vc- 

From the equations derived above, we can derive many important relationships. For example, 

bmce are independent of r, we can immediately conclude that the ratio 11 /n+ is also independent of r. This is 
an important conclusion from our calculation. From the above equations we can also derive 

K- = r'/'^-y-j. (70) 



(/-)4/3 
2 



— = r r- 

K+ \I 

The total energy cascade rate can be written in terms of E^{k) as 



(71) 



n = i(n+ + n-) = ^(/+ + r){E+{k)f'^k'^/^ (72) 

Since is independent of r, 11 is a linear function of r. When we apply the above formula to the solar wind stream 
with r = 0.07, we find that K+ = 2.12 and R- = 0.85. 

As mentioned in the previous section, the observed temperature evolution was studied by Verma et al. jl^ and Tu 
[2^. For streams with CTc — ^ 1, Verma et al. [2^ had assumed that K'^ = K~ = K independent of CTc, and derived 
the total turbulent dissipation rate to be 

n = ^(i?+(fc))^/^fc^/^. (73) 

Clearly the assumption that = K^, as well as the above formula H73|l is incorrect. Hence, the calculation of 
Verma et al. needs to be modified. The substitution of the parameters r and in our formula (|72l) gives us 
an estimate of the turbulent heating that is an order of magnitude higher than the observed overall heating in the 
solar wind [2^. Some of the resolutions of this paradox are: (1) the assumption that the solar wind has reached 
steady-state is incorrect, and the formula H72|l is inapplicable to the solar wind streams with large or (2) the 
constants calculated above will be modified significantly by the mean magnetic field, anisotropy, and helicity etc. 
In case of the former, one needs to understand the nonequilibrium evolution of MHD turbulence, while in case of the 
latter, the field theoretic calculation has to generalized in presence of mean magnetic field and helicity. Both these 
generalizations are beyond the scope of this paper. 



IV. SUMMARY AND CONCLUSIONS 



In this paper we have theoretically calculated various energy cascade rates in the inertial range of nonhelical MHD 
turbulence. Our procedure is based on field-theoretic approach. Using the steady-state condition we also calculate the 
energy supply rate from the large-scale velocity field to the large-scale magnetic field. For simplicity of the calculation, 
we have taken two special cases: (1) ac = 0; (2) CTc 1. Throughout the calculation we assume that the velocity 
modes at large length-scales are forced. 

We will first summarize the results for CTc = case in d = 3. The cascade rates HJ^> , n|j^, , H^>are approxi- 
mately the same for ta in the range of 0.5 — 1, but the flux HJJ^is rather small. The sign of H^^is positive, indicating 
that the ME cascades forward, that is from large length-scales to small length-scales. The large-scale magnetic field is 
maintained by the H^^ flux. We exploit this result to construct a dynamo model for galaxy. This result is discussed 
in paper II. 

Recently Cho and Vishniac (CV) performed numerical simulation of nonhelical MHD turbulence and arrived at 
the following conclusion based on their numerical results. In our language, their results for large ta can be rephrased 
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as (1) Il'^< « (2) n^^< ^j,^) « UB^; (3) n^< ^ {U - cB)B^, where U and B are the large-scale velocity and 
magnetic field respectively, and c is a constant. When we compare our theoretical findings with CV's result, we 
find our results can explain CV's first and second results, but they are only partly consistent with the third result. 
From Eq. H44|l it can be easily seen that n,"> depends on the KE in the same manner as in fluid turbulence. Hence, 
njl> ~ U^/L, a result consistent with the first result of CV. Using n"^^ ^^^^ — + + and the definitions 
of F's [Eqs. ()46I47|I ] we can easily show that 

+ (74) 

where ... represents a constant. We estimate the above equation in the large limit (i?" ^ E''). In this limit, 
n w U^/L. Hence, to a leading order in r^^ 

%^<+b» « UB^L (75) 



From Eqs. (|44I47|) . we also conclude that 



U B 
'UL 



n«< « {...r-/ + ...r-/)Ii « ( - ) B^ (76) 



Note that the first part of the above equation matches with CV's first part, but the second part of HJ^^ differs from 
CV's result by a factor of B/U . Since B/U « 1 at steady state, it is difficult to differentiate our results with that of 
CV. On the whole, our theoretical calculation is able to explain the numerical results of CV. 

For d = 'd the Kolmogorov's constant K is approximately constant and is close to 1.5 for all ta greater that 1/2, 
same as that for fluid turbulence {rA = oo). This result implies that the total cascade rate does not change appreciably 
under the variation of ta (since H oc K'^/"^). The cascade rates vanish near ta = 0.25; this result is in the expected 
lines because MHD equations become linear in = limit. Comparison with the past results shows that our result 
differs from that of Verma and Bhattacharjee's calculation [sj where Kolmogorov's constant changes significantly 
with the variation of ta- Note, however, that our procedure described here is an improvement over that of Verma 
and Bhattacharjee, where they had assumed a wavenumber cutoff for the self-energy integral for curing the infrared 
divergence problem. They had also assumed a specific type of self-energy matrix which can be shown to be correct 
only in some regime. 

When we vary d, we find that for large d, H^:< = n^< = n|;< = H^^ . In addition we also observe that Kolmogorov's 
constant MHD turbulence increases with dimensions as d^/^ . The same variation is observed for fiuid turbulence |24| . 
This result indicates that the cascade rates decrease in higher dimensions. We could calculate fluxes for d > 2.2 
because the RG fixed point is unstable for dimensions lower than 2.2 (Verma However, the RG fixed point for 

fluid turbulence is stable for c? = 2, and the Kolmogorov's constant in the inverse cascade regime of 2D fluid turbulence 
comes out to be 6.3. For this computation, the renormalized viscosity v* was taken as -0.6. It is interesting to note 
that Dar et al. find negative KE flux (H">) in their 2D MHD turbulence simulation; this is reminiscent of 2D 
fluid turbulence. 

In the other extreme limit Uc ^ 1 and r^i = 1, we find that Kolmogorov's constants and are not equal, 
and the ratio K~ /K^ depends very sensitively on r = E~ (k) / E'^ {k) . Both the fluxes H^, and also the total flux H, 
are proportional to r. The flux ratio H^/H+ is found to be independent of r. We also discuss the implications of our 
flux results to the heating of the solar wind. 

In this paper we have restricted ourselves to nonhelical turbulence. Helical MHD turbulence is very important 
specially for the growth of magnetic energy (dynamo). The energy fluxes for helical MHD have been discussed in 
paper II. The study of the effects of mean magnetic field using field theory has been relegated for future. 
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APPENDIX A: VALUES OF Ti 

The algebraic expressions for Ti{k,p, q) are given below. 

Ti(/c,p,g) = hP+^{k)P,aip)P^biq) 

kp{{d-3)z + {d-2)xy + 2z^ + 2xyz^ + x^z) (Al) 

T3ik,p,q) - hPr^kik)Pjaip)Pb{q) 

^ -k^{{d-2){l-y^) + z^ + xyz) (A2) 

T5ik,p,q) = -hP+Jp)P,aik)P,i,iq) 

= -kp{{d-3)z + {d-2)xy + 2z^ + 2xyz^ + y'^z) (A3) 

T7ik,p,q) = -hPjabiP)Pia{k)Pb{q) 

= -kp ((2 - d)xy + (1 - d)z + y^z) (A4) 

nik,p,q) = -hP+^{q)P,a{k)P,b{p) 

= —kq[xz—2xy'^z — yz'^^ (-A-5) 

Tn(fc,p,g) = -hPr^,{q)P,a{k)P,b{p) 

= —kqz{x + yz) (A-6) 

T2n{k,p,q) = -T2n-i{k,p,q) for n = 1..6 (A7) 

Ti3,i5ik,p,q) = kiMjab{k')Pja{p)Pib{q) 

= —kpyz{y + xz) (A8) 

Ti4:,i6{k,p,q) = kiMjab{k')Pjb{p)Pia{q) 

= k\l-y^)id-2 + z^) (A9) 

Ti7,i9ik,p,q) = kiMjab{p)Pja{k)Pibiq) 

= kpxz{x + yz) (AlO) 

Ti8,2Q{k,p,q) = -Ti4{k,p,q) (All) 

T'21,23(fc,P, ?) = kiMiab{q)Pja{k)Pjb{p) 

= -kpxy{l-z^) (A12) 
T22M{k,p,q) = -Ti2.{k,p,q) (A13) 

where k = p + q, and x,y, z are defined by 

p • q = ~pqx; q • k = qfcy; p • k = pfcz. (A-14) 
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Figure Captions 



Fig. 1. Various energy cascade rates of MHD turbulence. The illustrated wavenumber spheres contain u < and 
b < modes, while u > and b > are modes outside these spheres. The velocity fields is forced at large-scale. 
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TABLE I: The flux ratios computed by Dar et al. {d = 2, « 0, rA « 0.5) 

n^^/n ut</n ui</u nl</u n;^</n 

-0.13 0.68 -0.09 0.47 0.37 w4 



TABLE II: The computed values of energy cascade rates of MHD turbulence for various rA when d = 3 and cTc = 0. 



11 \ rA 


OUUU 


1 nn 
iUU 





-1 
1 


U.O 


U.O 


Trend 




1 


0.97 


0.60 


0.12 


0.037 


0.011 


decreases 


n;^</n 


3.5E-4 


1.7E-2 


0.25 


0.40 


0.33 


0.36 


increases 
then saturates 




-l.lE-4 


-5.1E-3 


-0.05 


0.12 


0.33 


0.42 


increases 
then saturates 


n^^/n 


2.7E-4 


1.3E-2 


0.20 


0.35 


0.30 


0.21 


increases 
then dips 


n«</n 


1.7E-4 


8.1E-3 


0.15 


0.47 


0.63 


0.63 


increases 
then saturates 


K+ 


1.53 


1.51 


1.55 


1.50 


1.65 


3.26 


approx. const 
till rA w 0.5 




1.53 


1.50 


1.29 


0.75 


0.55 


0.75 


decreases 



TABLE IV: The computed values Kolmogorov's constants for cTe — > 1 and ta = 1 limit for various r = E /E'^ {d = 2, 3) 



d r K+ I£_ 

0.17 1.4 1.4 

0.10 2.1 1.2 

3 0.07 2.7 1.07 

10"* 45 0.26 

10^^^ 4528 0.026 

0.1 1.2 2.4 

0.07 1.5 2.2 

2 0.047 1.9 1.9 

10"'^ 25 0.52 

10"^ 2480 0.052 



TABLE III: The 
rA = 1. 


computed values of energy 


cascade rates of MHD turbulence for 


various space 


dimensions d when Uc 


= and 


n\d 


2.1 


2.2 


2.5 


3 


4 


10 


100 


n^^/n 
u]:</u 




0.02 
0.61 
-0.027 
0.40 
0.37 


0.068 
0.49 
0.048 
0.39 
0.4 


0.12 
0.40 
0.12 
0.35 
0.47 


0.17 
0.34 
0.18 
0.31 
0.49 


0.23 
0.27 
0.23 
0.27 
0.50 


0.25 
0.25 
0.25 
0.25 
0.50 


K+ 




1.4 
0.69 


1.4 
0.72 


1.50 
0.75 


1.57 
0.79 


1.90 
0.95 


3.46 
1.73 
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FIG. 1: Various energy cascade rates of MHD turbulence. The illustrated wavenumber spheres contain u < and b < modes, 
while u > and b > are modes outside these spheres. The velocity fields is forced at large-scale. 



